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Abstract 

The annihilation of dark matter particles captured by the Sun can lead to a neutrino 
flux observable in neutrino detectors. Considering the fact that these dark matter 
particles are non-relativistic, if a pair of dark matter annihilates to a neutrino pair, the 
spectrum of neutrinos will be monochromatic. We show that in this case, even after 
averaging over production point inside the Sun, the oscillatory terms of the oscillation 
probability do not average to zero. This leads to interesting observable features in 
the annual variation of the number of muon track events. We show that smearing of 
the spectrum due to thermal distribution of dark matter inside the Sun is too small 
to wash out this variation. We point out the possibility of studying the initial flavor 
composition of neutrinos produced by the annihilation of dark matter particles via 
measuring the annual variation of the number of /x-track events in neutrino telescopes. 
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1 Introduction 



Growing evidence from a wide range of cosmological and astrophysical observations shows 
that about 82 % of the matter content of the universe is composed of Dark Matter (DM) 
whose exact identity is yet unknown {i.e., Pdm / {pdm + Pbaryon) = 82 % [Ij). In the literature 
[2] , various candidates for DM have been suggested among which Weakly Interacting Massive 
Particles (WIMPs) are arguably the most popular class of candidates [3j. The WIMPs are 
expected to propagate in the space between the stars and planets just like an asteroid which 
is subject to gravitational force from the various astrophysical bodies. In time a considerable 
number of WIMPs will interact with the nuclei inside the Sun and lose energy. If the velocity 
drops below the escape velocity, the particle will be captured by the gravitational potential of 
the Sun The DM particles will eventually thermalize inside the Sun. Having a relatively 
large density inside the core, the DM particles will annihilate with each other. Depending 
on the annihilation modes, two classes of the WIMP models can be identified: (1) DM pair 
directly annihilate to neutrino pairs (that is ui? or uu + iJi?). (2) DM pair particles annihilate 
to various particles whose subsequent decay produce neutrinos alongside other particles. 
In both cases, the neutrinos can be detected in the neutrino telescopes such as IceCube, 
provided that the cross section of the WIMPs with nuclei is of the order of ~ 10~^ pb 
or larger. Indirect detection of DM through registering the neutrinos has been thoroughly 
studied in the literature [5]. 

The oscillation length of the neutrinos due to 12- mixing can be estimated as 

^°"=A^~'"1° ""(l50fevj( Ami, )• 
which is of order of few percent of the distance between the Sun and the Earth. If the 
energy resolution of the detector {6E/E) is worse than 1 % and the width of the spectrum is 
larger than 6E/E, averaging the oscillatory terms is justified. However, for a monochromatic 
spectrum, dropping the oscillatory terms may lead to an error. Investigating the observable 
effects of these terms and information that they carry is the subject of the present letter. We 
shall demonstrate this effect by diagrams based on explicit calculation. We also demonstrate 
that integrating over production point does not justify dropping the oscillatory terms. 

The DM particles captured inside the Sun are non-relativistic and have an average velocity 
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of {3TQ/m£)MY^'^ — 60 km/sec. As a result, in the case that the DM pairs directly annihilate 
to neutrino pairs, the spectrum of neutrinos will be monochromatic. Considering the fact 
that Lose is of the order of the variation of Earth and Sun distance over a year, we expect 
these oscillatory terms lead to a significant variation of the number of events during a year. 
This expectation is similar to prediction of seasonal variation of beryllium line in the seminal 
paper by Gribov and Pontecorvo [B]. 

In the papers of reference [7], dropping the oscillatory terms, the oscillation of the neu- 
trinos on the way to the detectors has been studied. In papers [H HI [10] propagation of the 
monochromatic neutrinos has been numerically studied within the 3-z/ oscillation scenario; 
but, no emphasis has been put on the potential effects of the oscillatory terms. In |10] the 
potential effects of the oscillatory terms and the seasonal variations have been pointed out 
but since the emphasis was on the flavor blind initial composition, the effects had not been 
explored. 

In this paper we focus on the effects of oscillatory terms on the flux of monochromatic 
neutrinos from the Sun without assuming a democratic initial flavor content. Monochromatic 
flavor dependent neutrino flux can emerge within various models. A particular scenario 
leading to monochromatic neutrino flux has been recently worked out in [11]. The scenario 
is based on an effective coupling of form gaN(j)Va where is a Majorana neutrino and cf) is 
the scalar playing the role of the dark matter. Within this scenario, the main annihilation 
mode is DM + DM — )■ z/q, + 1/^3, z/^ + z/^ and the flavor composition can be determined by the 
flavor structure of the coupling g^. The scenario can be embedded within various models 
that respect the electroweak symmetry |12]. Another example of a model that can give 
rise to monochromatic neutrino spectrum with non-democratic flavor ratio can be found in 
|13j . We evaluate the magnitude of seasonal variation and discuss the condition under which 
such an amount of variation can be in practice established. We propose using the seasonal 
variation measurement as a tool to probe the physics of DM annihilation. 

In section |2l taking into account the matter effects inside the Sun, we derive the minimum 
value of the width of the spectrum leading the oscillatory terms to average to zero. We also 
show that uncertainty in the production point will not lead to vanishing of the oscillatory 
terms. In section |3l we investigate various sources that lead to widening of the spectrum of 
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monochromatic neutrinos. We then show that this widening is too small to make dropping 
of the oscillatory terms justifiable. In section |U we demonstrate the observable effects of the 
oscillatory terms in seasonal variation of the flux. In section El we anticipate the conclusions 
that can be reached based on different observational outcomes. The results are summarized 
in section |6l 

2 Effects of the oscillating terms 

In case that the DM pair annihilates to a neutrino pair and the physics governing dark matter 
annihilation is flavor blind, neutrinos of all three flavors can be produced in equal numbers. 
In this case, the neutrino oscillation cannot alter the flavor composition: if = = F^, 
— ^p- However, DM interaction can in general be flavor sensitive so neutrinos of 
different flavors can be produced with different amounts. For example, it has been recently 
suggested that there might be a correlation between flavor structure of DM couplings and the 
neutrino mass matrix [11]. In fact, it is possible that a coherent mixture of different flavors 
are produced at the source which means that the density matrix can be non-diagonal in the 
flavor basis. The density matrix, being a Hermitian matrix, can always be diagonalized. Let 
us denote the eigenstates of the density matrix at the source with \va) which are in general 
linear compositions of \ve), and |z/^). 

Let us consider such a state with momentum p produced inside the center of the Sun. 
This state after propagating a distance L evolves into 

|z/a;p; L) = aai{L)\l]p) + aa2{L)\2;p) + aa3(L)|3;p) (2) 

where \i;p) is the neutrino mass eigenstate (in the vacuum). Let us define (Ap*-');^^ as the 
minimum value of the energy interval for which 

r [afsi (0) {ai3j (0) )* (aa^ (L/ c) ) *aaj (L/ c)] |p+(Ap'^ )„, 
L [ai3i{0){apj{0))*{aai{L/c))*aaj{L/c)] \p 

where a P and L is the distance between the Sun and the Earth. Suppose that the 
spectrum of the neutrino fiux is almost monochromatic with width Ap ^ p. If the width of 
the spectrum is smaller than {Ap'^^)iim, the average of the oscillatory terms will in general 



= 27r , (3) 
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Figure 1: The dependence of {Ap^'^)iim/p and {Ap^^)iim/p defined in Eq. ([3]) on tlie momentum 
of neutrinos. The behavior shown in these graphs does not change by varying ^13 in the allowed 
range ^13 < 10 ° or changing the flavor composition. 

be nonzero. In this case, if the energy resolution of detector, 6E, is larger than Ap, the 
corresponding oscillatory terms will lead to an error in deriving the initial flux. Let us 
evaluate the numerical value of {Ap^^)iim- In vacuum, {Ap^^)iim can be estimated as 



(Ap 



lim 



P 



2-71 p 



(4) 



so taking L = 1.5 x lO^'^ cm we find {Ap^^)iim/p = 0.01 (p/100 GeV) and {Ap^^)iim/p = 
0.0005 (p/100 GeV). Notice that this result is independent of a and (3. Of course, the 
neutrinos produced in the center of the Sun have to pass through the matter inside the Sun 



before reaching the detectors. Figs. ([U-a) and ([U-b) show {Ap )iim/p and {Ap^)iim/p versus 
the neutrino momentum, taking into account the matter effects for various values of 613. 
To draw these figures, we have numerically obtained the evolution of the neutrino states of 
different flavors. Considering that the density of matter is sizeable only inside the Sun with 
radius Rq ~ 7 x 10^° cm ~ Lose, we expect the numerical values of {Ap)iim/p in the presence 
of matter to be almost similar to that in vacuum. In fact, the results shown in Fig. [1] fulfill 
this expectation: The order of magnitude of {Ap^^)iim/p is the same as that in vacuum and 
moreover independence of {Ap^^)iim/p from flavor {i.e., a and fi) is maintained. 

Similarly to Eq. ([2]), we can write 

|z/„;p;L) = ac,i{L)\l;p) + a„2(/^)|2;p) + 0^,3(^)13; p), (5) 
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where is the antineutrino mass eigenstate (in the vacuum). Due to matter effects 

inside the Sun, even in the absence of the CP- violating phase 6, aai{L) ^ aai{L). However 
since averaging of the oscillatory terms is mainly due to the propagation of neutrinos and 
antineutrinos in the large empty space between the Sun surface and the Earth, the value 
of {Ap''^)iim/p for neutrinos and antineutrinos will be similar. Direct numerical calculations 
confirm this claim. 

Another effect that may lead to averaging of the oscillatory terms is the difference in the 
production point; i.e., the difference in baseline. The dark matter particles are distributed 
in a radius of 



where T and p are respectively the temperature and density inside the Sun center. 
is the Newton gravitational constant. Inserting the numerical values, we obtain r£,M = 
2 X 10^ cm(100 GeV/mDAiy^^ = O.OO3i?0(lOO GeV/m^M)^/^- Inside the volume around the 
Sun center with r < tdm, the density is high such that Ve = y^GpNe ^ t^m^^lv^ so z/g and 
Vf. in practice correspond to energy eigenstates. As a result, i/J within this volume does not 
go through oscillation. On the other hand, the oscillation length corresponding to the 2-3 
splitting, 47rp/Am^^^ ~ 10^° cm is much larger than r^u so averaging out the corresponding 
oscillatory terms is not justified. Let us define P^/s as the average of the oscillation probability 
over production point with keeping the oscillatory terms. Moreover, let us define (-Pa^) as the 
same quantity with dropping the oscillatory terms. If the ratio Pap/iPap) equals to one, it 
means when we integrate over the production point, the oscillatory terms average to zero. We 
have numerically calculated Pap/ {Pap) and displayed it in Fig. O As seen in Fig. [2], Pe^/ {Pefi) 
and Pfifj,/ (Pfifj) can substantially deviate from one which means dropping the oscillatory terms 
for monochromatic spectrum is not justified. Thus, the numerical calculation demonstrated 
in Fig. [2] confirms our simple analysis made above. Another interesting point is the significant 
sensitivity to ^^13. This is understandable because for ^13 > 0, the 1-3 resonance in the Sun 
can play a role. We repeated the same numerical analysis for the case that the initial neutrino 
state is a coherent combination of the neutrino flavor states and reached the same conclusion. 




(6) 
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(a) (b) 
Figure 2: The dependence of the ratios Pei_i/{Pe^) and P^^/{P^^) on the energy of the produced 
neutrinos in the annihilation of the DM particles inside the Sun. P^p is the oscillation probability 
Ua — ?• 1^/3 obtained by integrating over the production point of neutrinos inside the Sun taking into 
account the oscillatory terms. To perform this analysis, we have taken the neutrino mass scheme to 
be normal hierarchical. The (Pa/s) is the same quantity without taking into account the oscillatory 
terms. Deviation of the ratio Pai3/{Pai3) from one is a measure of the significance of the oscillatory 
terms. 

3 Width of the spectrum of neutrinos 

The spectrum of neutrinos from direct annihilation of dark matter particles at rest will be 
monochromatic. In reality, however several thermal and quantum mechanical effects lead to 
widening of the spectrum. Let us consider them one by one. 

Thermal widening: The dark matter particles can fall in gravitational potential of the 
Sun and despite the very low cross section can eventually be thermalized. The temperature 
in the Sun center is about 1 keV so the average velocity is given hj v = (3T/m£)MY^'^ — 
60 km/sec(100 GeV/m^iAf)^/^. The widening due to this velocity is 

^ ^ ^ ~ 10- (-^] (7) 

E c Vl-3 keVy V ^DM J ' ^ ' 

Decoherence due to gravitational acceleration: The DM particles become accelerated in 
the gravitational field of the Sun so neutrinos from the annihilation of a DM pair will have a 
coherence length smaller than they would if the acceleration did not exist. We follow the same 
logic as in [2] to calculate the wavepacket size. Notice that this is a quantum mechanical 
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:y'^5xio-"fHI^V^ (8) 



effect originating from the uncertainty principle. Let us hypothetically divide the path of 
the DM particles to successive segments of size vAt. At each segment, the wavepacket of the 
emitted neutrinos will have a width of Ap ~ 1/Ar. During At, the velocity of DM particles 
is changed by ATcGprAr/S so the average momentum of the neutrino wavepackets emitted 
during the successive segments will differ by 6p ~ AirGprArmDM /<i- Now, if 6p ^ 1/Ar, 
the two successive wavepackets will form a single wavepacket with length 2Ar. However, for 
6p > 1/Ar, the two successive wavepackets are incoherent. That is the coherence length can 
be found by equating Ar = l/6p which leads to 

Ap _ UnGprpM Y^'^ ^ m-u /^lOO GeVV^^ 
p V 3mz)M 

Thus, widening of the spectrum due to acceleration is quite negligible. The width of the 
spectrum of the neutrino flux is dominated by thermal fluctuations rather than quantum 
mechanical widening. 

Natural width of neutrino wavepacket: The annihilation time sets a natural limit on the 
wavepacket size of the produced neutrinos. That is 

p niDM 

where udm is the density of dark matter in the Sun center, (aannv) is determined by the DM 
abundance in the universe: {aannv) ~ 10~^^ cm^. Evaluating udm is more model-dependent. 
Udm cannot be larger than Cr0/(47rr|,^^/3) where C is the capture rate of DM particles 
by the Sun and r© is the Sun lifetime. r£,M determines the size of volume around the Sun 
center where the DM density is relatively high (see, Eq. ([6])). Annihilation of course reduces 
the DM number density but let us take udm ~ Ctq/ {Aixr^-^f^J?)) to obtain a conservative 
estimate for the natural width. The capture rate is given by 



^ PDM ( Mq . 



where pdm = 0.39 GeV cm~^ [16j and v^m ~ 270 km sec~^ [T7] are respectively the local 
density and velocity of DM particles in our galaxy; Mq and are respectively the Sun and 
proton masses. Of course, o'DM-nudeon is unknown but if the interaction is spin-dependent, 
it can be as high as C(pb) [TS]. The maximal possible capture rate is therefore O[10^^ sec^^]. 
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Inserting the numerical values, we find that 

nDAlivCann) in-38 

~ iU 

so even with overestimating nDM{v<Jann), the natural lower bound is too weak and Ap/p will 
be dominated by thermal widening; i.e., Eq. ([7]). 

Widening due to scattering: At energies higher than 100 GeV, some of the produced 
neutrinos can undergo scattering before leaving the Sun. Neutrinos undergoing charged 
current interactions convert into charged leptons which are absorbed in the matter and do 
not contribute to the neutrino flux. An exception is of course z/,- — )■ r because the subsequent 
decay of the tau lepton regenerates high energy Ur. On the other hand, neutrinos undergoing 
neutral current interactions are converted to another neutrino with somewhat lower energy. 
Thus, the neutrino spectrum emerging from the Sun surface is composed of a sharp line 
at Eu = moM superimposed over a continuous spectrum with Ey < muM- The ratio of 
neutrinos with E^ c:^ m^M '^o those with E^ < moM depends on the neutral current cross 
section which itself depends on the energy of neutrinos before scattering. This energy is 
in turn determined by moM- The mean free path of the neutral current interaction of the 
neutrinos with energy 100 GeV in the Sun center is 

5 X 10^5 cm-3\ /1.3 X 10^3^ cm^X 
^0 /V ^NC ) 

where we have used the data from [19]. Considering the fact that the matter density in 
the Sun falls with radius as e~^^'^^'^^'^\ the ratio of neutrinos undergoing neutral current 
interactions should be of the order of Q.IRq/ Inc — 5%. The ratio increases with energy 
and at E^ = 500 GeV reaches 35 %. We restrict our analysis to the case m^M < 500 GeV 
which is also theoretically and phenomenologically motivated. In this range, the sharp line 
in the spectrum is quite pronounced so we will consider only this sharp line and will neglect 
the softened part of the spectrum. If m/jM turns to be greater than about 500 GeV, this 
analysis should be reconsidered taking into account the softening effects due to the neutral 
current scattering. 



iNc = — - — = 1-5 X 10^ km 

no(TNC 
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4 Observable effects of the oscillatory terms: seasonal 
variation 



For a monochromatic neutrino flux, the rate of yU-track events in IceCube as a function of 
time can be estimated as 



dt J [L(t)]2 

h + 

where integration is over the volume inside the Sun where neutrinos are produced. Fj^^ 
(Fp^) is the flux of Ua (j^a) produced in unit volume, pp and p„ are respectively the number 
densities of protons and neutrons in the ice. cr?'?, cr^*^, o"?^ and a?*! are the cross-sections 
of the charged current interactions of and z/^ with proton and neutron, respectively, w and 
w are suppressions factors respectively due to the absorption of neutrino and anti-neutrino 
fluxes in the Sun. For the energies that we are interested, E„ > 100 GeV, w and w do not 
depend on the flavor but in general, w ^ w. is the muon range in the detector which is 
the same for muon and anti-muon [2(J|. Af.ff{9\t\) is the effective area of the detector which 
depends on the angle between the neutrino momentum and the axis of the array of PMTs in 
detectors, 9. Because of the tilt of the rotation axis of the Earth, this angle changes as the 
Earth moves in its orbit around the Sun. L{t) is the distance between the Sun and the Earth 
which varies about 3 % during a yea. Fiually, P„„ aud P,, are respectively tite oscillatiou 
probability of and Va ^ ^^u These probabilities can be numerically derived from 

the evolution of neutrino states: 

= [ — 1^ — + diag(V;, 0, Q)\^^\va) (11) 



and 



I- 



diag(K,0,0)]^,|z><,) (12) 



dt ' 2p 
where K = V^GpN^- 

•^As discussed in the beginning of section [2l Va can be a coherent combination of different neutrino 
flavor eigenstates that diagonahze the density matrix. We perform our numerical analysis for the case that 
Ua corresponds to a given flavor. However, as discussed, our results apply to the general case, too. 
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The following remarks are in order: 



• In Eq. f lTU]) . we have neglected the subdominant contribution from z/,- — )■ r — )■ which 
is suppressed by Br(r — )• ^uu) = 17 %. 

• In the case that the spectrum is continuous, the total flux has to be replaced with 
differential flux and an integration over energy has to be taken. 

• Due to the tilt of the Earth rotation axis, during Autumn and Winter in southern 
hemisphere, neutrinos entering the IceCube pass through the mantle of the Earth 
before reaching the detectors. However, we can neglect the oscillation of the neutrino 
inside the Earth mantle. The reason is that Am'^j/p ^ Ve so the effective mixing inside 
the Earth vanishes and no oscillation takes place in the constant density of the Earth. 
(Notice that although Ve in the Sun center where neutrinos originate is much larger 
than Ve inside the Earth mantle, Ve in the Sun surface is smaller than that in the Earth 
mantle. That is while the density in the Earth mantle is almost constant. Because of 
the difference in profile, the matter effects in the Sun and Earth are different.) 

• Production can be either lepton number conserving {i.e., DM + DM — )■ z/ + z/) or 
lepton number violating {i.e., DM + DM — )■ u + u or DM + DM — > z/ + z/). In the 
former case, the flux of neutrino and antineutrino will be obviously the same. In the 
latter long as the CP is conserved in annihilation, processes DM + DM — )■ z/ + z/ 
and DM + DM u + u take place with the same rate. Even if CP is violated, any 
asymmetry in the z/ and z/ fluxes will be a subdominant effect resulting from the 
interference of the tree level and loop level contributions. Thus, we can safely take the 
initial flux of neutrinos and antineutrinos to be the same. 

Let us define 

f''^^\dNjdt) dt 
N{to, At) = ^ . (13) 

j::^'''Aen{e[t])l[L{t)fdt 

Notice that if the oscillatory terms average to zero, iV(to, At) will be constant and indepen- 
dent of to and At. Let us therefore define 

^ 7Vfa.Af)-A^fa + AMyear-Af) 

N{to, At) + N{to + At, 1 year - At) 

11 



(GeV) 


A(20March, 186days) 


A(3April, 186days) 


^13 = 0° 


^13 = 10° 


^13 = 0° 


^^13 = 10° 


NH 


IH 


NH 


IH 


NH 


IH 


NH 


IH 


100 


18 % 


18 % 


9 % 


11 % 


12 % 


12 % 


6 % 


7% 


300 


57% 


57% 


37 % 


42 % 


60 % 


60 % 


39 % 


43 % 



Table 1: The values of A(to, At) defined in Eq. ( IT^ for both normal and inverted hierarchies 
of the neutrino's mass spectrum. The initial flux is taken to be composed of z/g and with 
two different values of energies Ei, = 100 and 300 GeV. On to = 3 April, the earth reaches 
the point on its orbit depicted in Fig. |3] and to = 20 March corresponds to the spring equinox. 



E, (GeV) 


A(20March, 186days) 


A (3 April, 186days) 


^13 = 0° 


^13 = 10° 


^13 = 0° 


^^13 = 10° 


NH 


IH 


NH 


IH 


NH 


IH 


NH 


IH 


100 


9 % 


6 % 


4 % 


1 % 


7% 


4 % 


3 % 


0.3 % 


300 


12 % 


7% 


6 % 


19 % 


13 % 


7% 


6 % 


20 % 



Table 2: The same as Table [T] except that the initial flux is taken to be composed of i/^ and 



In the absence of the oscillatory terms, A(to, At) vanishes. Deviation of A from zero is a 
measure of the strength of the oscillatory terms. For a continuous spectrum, A vanishes 
because the oscillatory terms average to zero so A determines if there is a sharp feature in 
the spectrum. 

Tables [1] and [2] show values of A for different values of mDM{= E^), 613, (to. At) and for 
both Normal Hierarchical (NH) and Inverted Hierarchical (IH) neutrino mass schemes. To 
fill the Table [H we have taken the initial fiux at source to be composed of merely z/g and 
Ue and in the case of Table [21 we have taken the initial fiux to be composed of only z/^ and 
Ufj,. For one set of the data, we have chosen (to = 20 of March, At = 186 days). Twentieth 
of March corresponds to spring equinojQ and to + At corresponds to autumn equinox. For 



* Spring equinox is called "Norooz" in Farsi and is the beginning of the year in the Iranian calendar. 
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Figure 3: Positions of the Earth on to = 3 April and tQ + At = 6 October, where At = 186 days. 

another set, we have chosen to = 3 of April and At = 186 days which correspond to the 
points in the orbit indicated in Fig. [3J 

We have taken two values of ^^13 (^^13 = and 6*13 = 10°) and set the CP-violating phase, 6, 
equal to zero. We have calculated A{to, At) for the case of exactly monochromatic spectrum 
as well as for a narrow Gaussian with width given by the thermal fluctuations of the dark 
matter [see, Eq. ([7])]. As expected from the previous section, the difference in A for two cases 
(monochromatic versus Gaussian with the width given in Eq. ([7])) is quite negligible. In both 
cases, A significantly deviates from zero which means seasonal variation due to oscillatory 
terms is significant and potentially measureable. 

Several features are obvious in the figures of Tables [1] and [2l As can be seen from Table [H 
in the case ^13 = 0, the values of A(to, At) (for both to = 20 March and to = 3 April) are equal 
for NH and IH. This equality is a consequence of the fact that for ^13 = 0, the contribution 
of Us to z/g is zero so both the NH and IH cases are equivalent. By comparing the figures 
in Table [1] with the figures of Table [2] it can be seen that the value of A(to, At) is typically 
larger for the annihilation of DM particle to electron neutrinos (DM+DM— VeVe)- Thus, if 
the measurements of neutrino telescopes (such as IceCube) show a large value for A(to, At), 
the annihilation of DM particles to z/g (and z/g) will be favored. 

In the real experiments such as IceCube, measurement of A(to, At) will be more tricky 
and statistical errors and background events have to be taken into account carefully. The 
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statistical error of the ratio A (to, At) in Eq. ( !T^ is given by the following formula 

2 N{to, At)N{to + At, lyear - At) 



where we have inserted 6N = V N. To evaluate 6 A, we should estimate what is the maximum 
number of events per year allowed within the present bounds. The bound from present 
observation on the total number of muon tracks depends on the shape of spectrum of the 
neutrinos. This is understandable because the detection threshold of AMANDA and its 
augmented version IceCube are not exactly the same, especially once DeepCore is added. In 
[2T] . an analysis has been made for the spectrum corresponding to annihilation into W'^W" 
as well as into r+r^ for various values of the DM mass. The result is that if the bound is 
saturated, depending on the DM mass or the shape of the spectrum, a km'^ scale detector 
can observe between 500 to a few 1000 events per year. The spectrum for annihilation into 
neutrino pair is harder so the bound will be stronger. Detailed analysis is beyond the scope of 
the present paper. Taking total number of 400 /year/km^ however sounds reasonable. With 
400 muon tracks, the statistical error, 6 A, from Eq. f|T5|) . is less than 0.05. From Eq. (|9]), 
we find that this corresponds to nucleon DM cross section of 10~^pb which is well below the 
bound on spin-dependent cross section. Typical values for A in Tables [1] and |2] can therefore 
be established after a few years of data taking {i.e., 6 A ^ A). In practice, the measurement 
will suffer from reducible and irreducible backgrounds. A well-known irreducible background 
comes from the "solar atmospheric neutrinos". That is neutrinos produced by interaction 
of cosmic rays with the atmosphere of the Sun which amounts to about 10 events per year 
[22] . The signal can be larger than this background by one order of magnitude. Moreover, 
the solar atmospheric flux, its energy spectrum as well as its flavor composition can be 
calculated so its effects can be subtracted. Two other sources of backgrounds exist for the 
/i-track events from the Sun: the atmospheric neutrinos and the atmospheric muons. From 
the spring to autumn equinoxes the Sun is below the horizon at the IceCube site and the 
neutrinos coming from the Sun produce up-ward going yU-tracks. The atmospheric muons 
will be absorbed in this period of time by the surrounding material of the detector. This 
results in huge suppression of atmospheric muon background such that it can be completely 
neglected. During the autumn to spring equinoxes, the Sun is above the horizon for IceCube 





(15) 



A^(to, At) + A^(to + At, lyear - At) 
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and the number of atmospheric muons increases drastically. However, by using the outer 
parts of the IceCube as a veto, the atmospheric muon background will be quite low at the 
DeepCore [23] (the surrounding area can veto atmospheric muon events up to one part in 
10^). Let us now discuss the muon-track events induced by atmospheric neutrinos. As 
shown in appendix C of Ref. [23j, thanks to the high angular resolution of the IceCube for 
yU-tracks, by focusing on a cone with half angle 1° around the direction of Sun, the number 
of background through-going /i-tracks induced by atmospheric neutrinos can be reduced to 
~ 6 yr^^ for the whole energy spectrum above the IceCube energy threshold of the yU-track 
detection (number of background contained //-track events in the DeepCore is ~ 2.5 yr^^). 
Thus, the background events from the atmospheric neutrinos are not also a limiting factor 
in the measurement of A(to, At). 

In the previous section, we observed that P^,^ is quite sensitive to ^13. The reason is that 
for = 0, no 1-3 resonance in z/g — )■ takes place but for ^^13 > 0, such a resonance plays a 
significant role. From the tables, we observe that A is also quite sensitive to ^13. Considering 
this high sensitivity, it is tempting to entertain the possibility of extracting ^13 from A but 
absence of knowledge on the initial flavor composition renders the method useless. 

5 Discussion of results 

Throughout this paper we have focused on a monochromatic spectrum [i.e., a very narrow 
line) of neutrinos from the direct annihilation of dark matter particles to neutrino pairs inside 
the Sun. In reality, as we discussed earlier, the neutral current interactions inside the Sun 
will smear a part of the sharp line into a continuous spectrum with lower energy. Moreover, 
along with DM + DM -^^v + ''v, DM particles can annihilate into rf and other particles 
whose subsequent decay will lead to a continuous spectrum with lower energy. Thus, the 
spectrum will be composed of a sharp line superimposed on the upper edge of a continuous 
background. Our discussion holds valid about the sharp line part of the spectrum. 

In neutrino telescopes such as IceCube, the direction of /i-track can be reconstructed by 
amazing precision of 1° [25] which means neutrinos from the Sun can be singled out. In 
practice, the measurement of the spectrum is going to be challenging especially that a large 
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fraction of muons are produced by interaction of outside the detector and lose energy 
before entering the detector. However, if the statistics is high enough, the spectrum can 
be reconstructed by measuring the energy of contained muons. Let us consider different 
situations one by one. 

In case that the spectrum is not reconstructed, the measurement of A gives invaluable 
information on the shape of the spectrum. From Tables [1] and [21 we see that A can be 
quite sizeable for both normal and inverted hierarchical mass schemes. Thus, if the number 
of events is a few hundreds, deviation of A from zero can be established. Since we have 
divided the number of events by j L^'^Aejfdt in the definition of A^(to, At), a deviation of 
A from zero indicates that the oscillatory terms in oscillation probability do not average 
to zero (see, Eqs. fll3|14p ). This in turn shows that there must be sharp features in the 
spectrum, originating from direct annihilation of DM pairs to neutrino pairs. To reach 
such a conclusion, the possibility of other seasonal modulation (like detector performance or 
seasonal variation of the background [2S]) has to be of course subtracted. 

In case that the spectrum is reconstructed, A again provides invaluable information. If 
the spectrum contains a sharp line, we in general expect A to be nonzero. If a sharp line is 
observed in the spectrum but A turns out to be zero, the most natural explanation is that 
the initial flavor ratio is F^^ = F^^ = F^^ which in turn means the physics of DM annihilation 
is flavor blind. This may be a unique way to learn about flavor composition as the cascade 
events most probably will lie below the detection threshold for moM ^ 500 GeV. In general, 
to analyze the properties of dark matter, seasonal variation. A, provides a powerful tool. 

6 Concluding remarks 

In case that DM pairs annihilate into neutrino pairs, there will be a sharp line in the spec- 
trum. We have shown that in the presence of such a line, the oscillatory terms in the 
oscillation probability do not average to zero and can therefore lead to a seasonal varia- 
tion of the number of events in neutrino detectors. We have shown that the main cause 
for widening of the line is the thermal velocity distribution of DM particles inside the Sun 
but the widening will be too small to lead to vanishing of the oscillatory effects. We have 
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demonstrated that the uncertainty in the production point cannot lead to vanishing of the 
oscillatory effects, either. 

We have defined an observable quantity, A, whose deviation from zero is a measure of 
the significance of the oscillatory terms {see, Eq. (fT^ ). We have shown that A can reach 
as high as 60 % so its deviation from zero can be established by a few hundred muon-track 
events. We have calculated the background and statistical error and have found that in 
the case that the flux is close to the present bounds, measuring A is doable. A contains 
invaluable information on the properties of DM particles. If A 7^ 0, even without performing 
the challenging energy spectrum reconstruction, we may conclude that there is a sharp line 
in the spectrum so the DM pairs have an annihilation mode into neutrino pairs. If the 
spectrum is reconstructed and a sharp line is identified but A turns out to be zero, a natural 
explanation is that at the source all three flavors are produced by equal amounts. This 
in turn means that DM annihilation is flavor blind. Considering that for low values of E^, 
{i.e., for rriDM 500 GeV), the cascade events will be below detection threshold of neutrino 
telescopes, A might indirectly provide a unique probe of flavor composition. 

As discussed, the width of the line is given by thermal fluctuations in the Sun center 
(r < Tom < O.OIRq). Even conventional solar neutrinos from thermonuclear processes are 
produced mainly in the outer layers so our information on this region depends merely on 
solar models. Our results are however robust against solar models as Ap/p scales as T^/^ 
{see, Eq. ([7])). In order for the oscillatory terms given by Am^g to be erased, T has to be 
10000 times larger than the value predicted by solar models which seems quite unlikely. 
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